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ABSTRACT 

Context. The quantitative study of the physical properties and chemical abundances of large samples of massive blue stars at differ- 
ent metallicities is a powerful tool to understand the nature and evolution of these objects. Their analysis beyond the Milky Way is 
challenging, nonetheless it is doable and the best way to investigate their behavior in different environments. Fulfilling this task in an 
objective way requires the implementation of automatic analysis techniques that can perform the analyses systematically, minimizing 
at the same time any possible bias. 

Aims. As part of the ARAUCARIA project we carry out the first quantitative spectroscopic analysis of a sample of 12 B-type su- 
pergiants in the galaxy NGC 55 at 1.94 Mpc away. By applying the methodology developed in this work, we derive their stellar 
parameters, chemical abundances and provide a characterization of the present-day metallicity of their host galaxy. 
Methods. Based on the characteristics of the stellar atmosphere/line formation code fastwind, we designed and created a grid of mod- 
els for the analysis of massive blue supergiant stars. Along with this new grid, we implemented a spectral analysis algorithm. Both 
tools were specially developed to perform fully consistent quantitative spectroscopic analyses of low spectral resolution of B-type 
supergiants in a fast and objective way. 

Results. We present the main characteristics of our fastwind model grid and perform a number of tests to investigate the reliability 
of our methodology. The automatic tool is applied afterward to a sample of 12 B-type supergiant stars in NGC 55, deriving the stellar 
parameters. Si , C , N , O and Mg abundances. The results indicate that our stars are part of a young population evolving towards a 
red supergiant phase. For half of the sample we find a remarkable agreement between spectroscopic and evolutionary masses, whilst 
for the rest larger discrepancies are present, but still within the uncertainties. The derived chemical composition hints to an average 
metallicity similar to the one of the Large Magellanic Cloud, with no indication of a spatial trend across the galaxy. 
Conclusions. The consistency between the observed spectra and our stellar models supports the reliability of our methodology. This 
objective and fast approach allows us to deal with large samples in an accurate and more statistical way. These are two key issues to 
achieve an unbiased characterization of the stars and their host galaxies. 

Key words. Stars: early-type - Stars: fundamental parameters - Galaxies: stellar content - Galaxies: individual: NGC 55 

1. Introduction ISimon-Dfaz & Stasinskall201 ll and iBresolin et al.ll2009l) . offer- 

ing at the same time access to chemical species that are not 



5^ , The latest generation of large telescopes has opened a wide accessible through HII region studies (like Si or Mg), and 
9; range of possibilities in the study of massive blue stars, al- ^t distances where nebular results strongly rely on techniques 
lowing for the first time analyses of resolved stars beyond the ^hjch need to be carefully calib rated (the so-called strong- 
Magellanic Clouds, even to nearby galaxies beyond the lim- ij^e methods, iPagel et all [Tgfl . Moreover, blue supergiant 
its of our Local Group. This new observational capability is ^tars present themselves as very promising distance indicators, 
especially important not only to reach a better knowledge of through the appli cation of the wind-mo mentum-luminosity re- 
the nature of these objects, but also to understand the cheim- lationship (WLR, iKudritzki et alJ ITgg?) and the flux-weighted 
cal_|nd dynamical evolution of their host galaxies (e.g. |U et alj gravity-luminosity relationship (FGLR, Kudritzki et al. 2003) . 
120091) . The last two decades in paiticular have witnessed the 

use of massive stars as reliable metallicity tracers, complement- !„ order to fully understand the nature of these objects, it is 
ing the classic approach based on H II regions (see for instance required to perform accurate analyses on large samples of mas- 
sive stars both in the Milky Way and external galaxies. Whilst 



Send ojfprint requests to: N. Castro (e-mail: norbertoOnoa . gr) current multi-object spectrographs are certai nly capable of pro - 

* Based on observations obtained at the ESO VLT Large Programme ducing such large collections of spectra (e.g. lEvans et al l l2005h . 

171.D-0004. the accurate modeling of their atmospheres is an intrinsically 
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complex task, involving non-local thermodynamic al equilibrium 
processes and strong stellar winds, hence requiring the use of 
highly sophisticated stellar model atmospheres. The important 
advances acc omplished in the mo delin g of massive blu e star 
atmospheres (ISantolava-Rev et al.|[T997l: iHillier & Milled[T99l 
iPuls et al]|2005h . together with the improvement in the compu- 
tational facilities provide us with the tools to overcome these is- 
sues. With regard to the analysis technique, several alternatives 
have been proposed by different authors to minimize the sub- 
jective component present in the widely used by-eye techniques, 
by introducing objective, auto matic and fast method s, such as; 
the methodology employed by Mokiem et al. ('2005') based on 
the genetic algorithm PIK AIA (Charbonneau 1995) . the grid- 
based method proposed bv iLefeveil (|2007|) or the principal com- 
ponen ts analysis (PCA) algorithm designed by lUrbaneia et alj 
(l2008h . In this work, we present a new automatic grid-based 
technique implemented over a grid of models calculated with the 
latest release of the model atmosphere/line formati on code fast- 
wind (ISantolava-Rev et al.l[T997t iPuls et alJl2005h . specifically 
designed and optimized for the study of O- and B-type stars in 
the optical and infrared range. The combination of the fastwind 
code and a grid-based technique enable us to perform, in an ob- 
jective and fast way, the analysis of B-type supergiants in the 
optical domain at low spectral resolution. This provides a very 
efficient way to handle the analysis of data collected by large 
spectroscopic surveys. 

The growing interest in the nature of massive stars and 
their host galaxies has propelled the development of several 



studies within different members of the Local Group, 
example: in the nearb y Magellanic C l ouds (iLennon et al.l 
RoUestonetaP 120031: iTrundle et all l2004t iKornetalJ 



Lennon et alj2005l:lTrundle & Lennoni2005llHunter et al.i2007l 
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2003t 



2005t 



ITrundle et all l2007h.M 31 (fVennetaLl 
l2002h. M 33 (iMonteyerde & HeiTerol 



19981; iMonteverde et alJ 120001: fUrbaneia et al.l l2005bi 



NGC 6822 (iMuschielok et al l 119991: IVennetalJ 
NGC 3109 (lEvansetalJ l2007h. W L M (IVenn et alJ 
iBresolin et al." '2006- 'Urbanei aet ap l2008l) or IC 
(Bresolin et al. 2007: Garcia et al.l l2010l) 



2011 



2001), 



2003t 
1613 
Moreover, these 
studies have exten ded to galaxies bey ond the Local Group , 
such as NGC 300 (IBresolin etal.lllooa lUrbaneia et al.l l2005at 
lKudritzkiet"ani2008h or NGC 3621 dBresolinet al.ll200ll) . the 
latter being at a distance of ~6.7 Mpc. 

The ARAUCARIA projecfl (P.I: W. Gieren, iGieren et al.l 

l2005l) is an ambitious project devoted to investigate the effects 
that the environment could have on different distance indica- 
tors. To that end, a number of nearby galaxies (NGC 6822, 
IC 1613, WLM, NGC 3109, NGC 55, NGC 247, NGC 300 and 
NGC 7793) have been targeted, both photometrically and spec- 
troscopically. An important part of this ESO long-term project 
is focused o n the young stellar p opulation of these gal axies (see 
for instance Evans et ani2007l or lUrbaneia et al.|[2008l) . and one 
of its main results has been the discovery and partial exploita- 
tion of the FGLR of B A supergiant stars as a distance indicator 
dKudritzki et al.ll2003h . Within the context of the ARAUCARIA 
project, we presente d the first quahtati ve analysis of massive 
blue stars in NGC 55 dCastro et al.ll2008[ hereafter C08) situated 
in the Sculptor fi lament at 1 .94 Mpc dPietrzvhski et al.l 120061: 
iGieren et al.ll2008k . We now present the first quantitative analy- 
sis on a sample of B-type supergiants in this galaxy. This is a key 
step for characterizing the prominent population of blue massive 



https://sites.google.com/site/araucariaproject/ 



stars in NGC 55, noted by IGieren et al.l (l2005l) . particularly in 
terms of their evolutionary status. 

The work presented in this paper is divided in two main parts. 
In the first one. Sect. |2] we present a detailed description of the 
tools we have designed for the analysis of massive blue stars 
at low spectral resolution. The goodness-of-fit criteria and the 
parameter space covered by our model grid, along with several 
tests to identify the limitations and the reliability of our method- 
ology are also presented. In the second part, we analyze 12 early 
B-type supergiant stars observed in NGC 55, using the previ- 
ously discussed tools (Sect. O. We will use the results to study 
the chemical distribution of the present day metallicity of this 
galaxy, in addition to constraining the evolutionary status of the 
analyzed massive blue stars, in Sect. |4] Finally, Section |5] pro- 
vides the final remarks and comments. 



2. A grid-based quantitative analysis technique 

The quantitative analysis of optical spectra of B- type supergiant 
stars i s based on well es tablished methods (e.g. ICrowther et al.l 
'200^ ISearleet ai1l2008l) . At high spectral resolution, the de- 
termination of temperature and surface gravity is based on the 
ionization balance of different ionization stages of the same el- 
ement (e.g. Si IV/Si III), and the fit to Balmer lines wings re- 
spectively (see lMcErlean et al]|1999l) . A slightly modified tech- 
nique is applied in the analysis of low spectral resolution data. 
Although the temperature and gravity criteria are the same, re- 
stricting the analysis to individual (metal) lines at low spec- 
tral resolution is unreliable. The best approach is to reproduce 
the main features present in the spectrum simu ltaneously, as it 
was suggested by lUrbaneia et al.l d2003L 12005 ab in the analysis 
of NGC 300 B-type supergiant stars. This technique has been 
also su ccessfully applied in the analysis of massive blue sta rs in 
WLM ("Bresolin etal."2006'), NGC 3109 dEvans et aLJIJOOl and 
IC 1613 (Bresolin et al. 2007). 

The complete spectral analysis consists of two steps. In the 
first one, the fundamental stellar and wind parameters are de- 
rived by using a fixed set of models. The determination of the 
chemical abundances is carried out in a second step by comput- 
ing tailored models. With the ultimate goal of performing quan- 
titative studies of low resolution spectra of OB-type supergiants 
at different metallicities in a systematic and objective way, we 
have implemented an automatic algorithm to determine the stel- 
lar parameters by identifying those models in our grid that mini- 
mize the differences with respect to the observed optical spectra 
(ax^ minimization). In the next sections we describe the main 
components of our automatic analysis method. 



2.1. FASTWIND grid of models 

The cornerstone of the analysis of massive blue stars is the grid 
of model atmospheres employed to reproduce the different fea- 
tures of the spectrum. Because of its high computational effi- 
ciency, w e have used the model atmosphere/line form ation code 
FASTWIND dSantolava-Rev et al.ll997l : lPuls et al.l2005l) . This code 
takes into account NLTE effects in spherical symmetry with an 
explicit treatment of the s tellar wind effects by cons idering a 
;6-like wind velocity law (ISchaerer & Schmutzl 119941) . and by 
ensuring a smooth transition between the pseudo-static photo- 
sphere and the inner wind layers. The main advantage with re- 
spect to other similar codes is the possibility of generating real- 
istic models in a short period of time, a crucial point for building 
large sets of synthetic spectra. 
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Fig. 1. Distribution of the models in the log Teff-log L/Lq plane 
(black dots). To illustrate the stars that can be analyzed with 
this gr id evolutionary tracks with rotation of iMeynet & Maeden 
(l2003h are plotted. The Humphreys-Dayidson Limit is also 
shown in the upper part of the diagram (solid line). 



2.1.1. Stellar parameters 

Each FASTWiND model is described by nine main parameters (we 
assume that the winds are homogeneous). Ideally, all of them 
should be considered free in our grid of models. However, this 
would imply the calculation of a very large number of models to 
explore the full parameter space. Therefore, we fixed and con- 
strained some of them based on the previous knowledge of the 
physics of these objects. This saves a significant amount of com- 
puting time, without introducing any relevant limitation in the 
analyses. A brief description of the criteria used to fix some of 
the parameters follows, as well as the specific range explored in 
each case. As a general rule, the boundaries of the model grid 
were defined aiming at avoiding observed stars too close to the 
grid's limits. 

- Effective temperature (Tetf). We computed models with 
temperatures between 9000 and 35000 K, in steps of 1000 K. 
This interval covers objects with spectral types ranging 
~A1 -08, for solar metallicity. 

- Surface gravity (log g). Whilst our main interest here is 
the analysis of supergiant stars, we extended the calculations 
to higher gravity values. In order to select t he gravities for 
each temperature, we considered the FGLR (iKudritzki et al.l 
l2003h . As shown by Kudritzki et al. (2008), this quantity is 
empirically related to the luminosity of normal blue super- 
giants 



M 



FGLR 



= (3.41 + 0.16) (log gF - 1.5) - (8.02 + 0.04) (1) 



where A^boi ^^ '■^^ bolometric magnitude, and 
log gf = log g - 4 1og(reff X lO""*). Supergiant stars 
of luminosity class la and lb, exhibit a range of log gf 
between ~ 1 .0 and 1 .5 dex. For objects with a temperature 
of 25000 K, this means log g ~ 2.6 - 3.1 dex. Our models 
were calculated covering a log gf range between 0.9 and 
2.5 dex, with an upper limit of log g = 3.7 dex (i.e., main 
sequence stars are not considered in this grid). 



Figure \T\ shows the location of our models in the 
Hertzsprung-Russell (HR) diagram, along with the evolu- 
tionary tracks from Mevnet & Maeder (2003j) for an initial 
equatorial rotation of 300 km s"' and solar metallicity. As 
can be deduced from this figure, our parameter selection cor- 
responds to stars with initial masses between ~ 8 - 60 Mq. 
Since we have constructed the grid with constant log gf val- 
ues instead of log g, our models define constant luminosity 
sequences. 

Radius (/?»). For each given pair [Tetf , log g] the radius was 
calculated from M/,o/ by means of the FGLR (see Eq. [TJ. 
Note that this relationship was observationally established 
for supergiant stars, but we have also used it for models that 
would represent stars that do not belong to this luminosity 
class. This, however, has no eff'ect on the analysis. 
Microturbulence (^). Three different values of the microtur- 
bulence, 7, 17 and 27 km s ', were considered for the calcu- 
lation of the model atmospheres. A larger number of values 
were used in the computation of the formal solutions. This 
procedure does not lead to inconsistencies in the spectrum 
as long as the value used in the formal solution does not de- 
part far from the value used in the model atmosphere. Hence, 
the formal solutions were calculated with ^ equal to 5, 7, 10, 
12 kmi ' (in the case of model atmospheres computed for 
7kms-i), 15, 17, 20, 22 km^"' (for Hkms"'), 25, 27 and 
30kms-i (for27kms-i). 

Helium abundance (He/H). The helium abundance by 
number is sampled with four points 0.05, 0.1 (solar), 0.15 
and 0.2. We note here that, although our lowest value is be- 
low the primordial He abundance and hence is physically not 
realistic, it was set to avoid boundary issues, as previously 
discussed. 

Metallicity (Z). Five values were used, from Z = 0.25 to 
1 .25 20 in step s of 0.2 5, with the solar references taken from 
lAsplund et al.l (l2009l) . For each metallicity, the elemental 
abundances (excluding He) are scaled by these values. 

Terminal velocity (Voo). For each model, this parameter was 
obtained from the escape velocity (Vesc) by using an em- 
pirical calibration based on the works by 'Kudritzki & Puis] 
(2000), Crowtheret al. (2006) a nd Markova & Puis (2 0oC 
According to th e studies of iLamers et al.l (1 1995b and 
IVink et all d 19991) there is a bimodal relationship between 
both quantities, with an abrupt change at the location 
of the so-calle d 'bi-s tability jump'. Contrary to this idea, 
ICrowther et al.l (l2006l) argued that there is no such break 
in this relationship. Rather, these authors propose a smooth 
transition around 20000 K. With the goal of emulating the 
empirical values, a monotonous trend was considered in the 
range of temperatures where the jump is located, 



n.io 

-^ = 1 11.65 log Teff - 47.62 
^''"^ 341 



Teff < 15 kK 

15kK<Teff <24kK (2) 

Teff >24kK 



To account for metallicity effects we assumed that the 
termi nal velocity scales with m etallicity as Voa(Z) 



Z°^^ (iLeitfaerer et al.lll992i see also IKudritzki & Pu"isll2000l 



IVink et al.ll200ll and lKudritzkill2002h 



Wind velocity law, j8. We adopted an empirical Unear re- 
lationship between ;6 and Teg, bas ed on results obtained by 
,Urbaneia ( 200J ) and lCrowtheret al. (2006) for Galactic B- 
type supergiants with temperatures in the range ~ 10000- 
~31000 K. Beyond those limits we used fixed values: 
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[3.60 reff<10kK 

j6 = |-1.40(reffl0-*) + 5 10kK<Teff <30kK (3) 
[o.70 reff>30kK 

Mass loss rate (M). For each combination of T^fi, log g, ^, Z 
and He/H, we conside red 3 different values for the mass loss 
rate. iPuls et al] (1 1996b showed that different combinations of 
mass loss rate, terminal velocity and radius can produce the 
same emergent synthetic profiles as long as the optical depth 
invariant, defined for smooth winds sls Q - M / (7?,Voo)' , 
remains constant. Since /?* and Voo are not free parameters 
in our grid, variations of M are equivalent to variations of 
Q. In order to use realistic mass-loss rates for each model, 
we applied the empirical relationship betwee n the wind mo- 
mentu m an d the stellar luminosity fo und by iKudritzki et al.l 
(Il995h and IKudritzki & Pulsl ( l2000h to define a proper Q- 
value. The wind momentum-luminosity relationship is de- 
fined by 



Table 1. Spectral features used in the determination of the fun- 
damental parameters. 



log D,nom = X log LJLq + D^ 

s 2.31 log LJLa + 15.94 



(4) 



where Dmom is the modified wind momentum (log D^om - 
log {MvooRJ )). The values for x and Do used in Eq.|4]were 
derived from a li near regression to data from previous studies 
(TCrowther et a ll 120061; iMokiem et all l2007bt iLefever et all 
l2007;,Markova & Pulsll20'08l) ^ 

For each set of stellar parameters three different Q values 
are considered: the one derived from Eq. log Q, and two 
others with log Q increased/decreased by 0.5 dex. Finally, 
we set a lower limit for M at 10"^ Mq yr"', since at this low 
mass loss rate the wind effects on the optical profiles are neg- 
ligible. 

The metallicity dependence of the mass-loss rate is ac- 
counted for with a power- law, M(Z) oc Z ' ". For th e exponent, 
we use the results from iMokiem et al.l (l2007bh that found 
m - 0.83 based on the analysis of Galactic, LMC and SMC 
stars. 



2.1.2. Atomic models 

The atomic models used in the calculations will play an im- 
portant role not only in the determination of stellar parame- 
ters but also in the computational time required per model. 
Detailed atomic models of HI, Hel-II dJokuthvl [200^, Sill- 
III- IV (N. Przybilla 2007, private communication) and of O II- 
III (lBecker&Butleilll988l: D. Kunze 1998, private communi- 
cation) are explicitly considered (see below) during the stel- 
lar parameters determination. The inclusion of O is critical 
since at low spectral resolution the line of Si IV 4116A and 
Sill 4128 - 30 A are blended with several Oil transitions. 
During the chemical abundance analysi s, detaile d models for 
CII-IIllE ber 1987: 'Eb er & Butleilll988h . N II-III (Butler 1984; 
iBecker & Butler 1989) and Mgll (K. Butler 1998, private com- 
munication) are also incorporated for the calculation of the tai- 
lored models. 



We note here that the other species are treated in an implicit 
way to account for blanketi ng/blocking efi'ect s. For further de- 
tails, the reader is referred to lPuls et al.l(l2005l) . 
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Notes. At low-resolution the lines of Sill 4128.07 A and Si IV 
4116.10 A could be blended with Oil transition. 



2.2. Stellar parameters determination 

To avoid a subjective and time-consuming procedure we have 
implemented a straightfo rward y^ technique, similar to the 
methodology proposed by iLefever et al.l (l2007l) . The observed 
spectrum is compared to a grid of synthetic models, in a number 
of relevant optical lines. To characterize the goodness-of-fit the 
differences are evaluated through the following the expression: 



"lines jT^ "v -^^ V Cr / 



(5) 



where «v is the number of wavelength points in the spectral 
line j, Vobs and yij are the observed and synthetic fluxes respec- 
tively (the index / runs over the set of models), and the un- 
certainty cr is estimated according to the signal-to-noise ratio 
(SNR). Eventually, the average of all the transitions is consid- 
ered, with all the lines having the same weight. The method was 
tested giving more weight to those lines that could have more 
impact on particular stellar parameters (e.g. silicon transitions in 
the effective temperature). The results revealed a better match 
on average when no extra weights were imposed. The selected 
lines are listed in Table[T] This selection is based on the observed 
range, the quality of the data and previous knowledge of model- 
ing these spectral features. 

The stellar parameters, and their uncertainties, are deter- 
mined through two steps, aiming at accounting for different 
sources of uncertainties: 

- From the x^ distribution generated by Eq. |5]we selected all 
those models with a x~ value below its x^ minimum plus 
15%. This percentage was not chosen arbitrarily, but us- 
ing a sample of high resolution spectra with SNR larger 
than 150 obtained as part of the IACOB|j spectroscopic 
database. We carried o ut their analyses by a classic method 
(iMcErlean et al.lll999l) and with our new x^ minimization. 
By comparing the results, we could identify the percentage 
that recovered similar errors in both methods. 
Those models with ?iX^ that fulfill this criterion were chosen 
to derive the stellar parameters. The values and their errors 
were calculated by averaging all these models, weighted by 
g-o V (normal distribution probability), and their standard 
deviations respectively. 



^ The lACOB spectroscopic database (Simon-Diaz et al., in prepara- 
tion) presently contains ~ 200 high quality spectra of O- and B-type 
Galactic stars {R = 46000, SNR> 150) 
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Fig. 2. Stellar parameters derived for the Galactic star HD 14818 (R - 1000 and SNR - 100). The panels show the parameters 
obtained in each Monte Carlo (MC) iteration, the average values are marked by a red solid line (see text for more details). The 
mean error obtained from the standard deviation in each simulation is marked by red dashed lines while the dispersion in the MC 
simulation is shown by blue dashed Unes. 



We investigated the effect that applying different percentage 
cuts would have on the derived parameters. For example, the 
difference between a 10% or 20% cut is very small when 
compared with the errors, with an increment of ~ 200 K in 
the effective temperature uncertainty. 

This method was complemented with a Monte Carlo simula- 
tion. Given the SNR of the spectrum, a random array of 100 
elements around the continuum position was defined. New 
re-rectifications of the stellar continuum were performed us- 
ing these shifts, then the stellar parameters and errors were 



re-calculated (through the step described before). This al- 
lows us to evaluate the impact that the SNR and the con- 
tinuum rectification have on the results. 



Both steps produce a set of errors in the final stellar parame- 
ters. We are aware that they are not independent and both estima- 
tions are linked to the spectral SNR and the intrinsic uncertain- 
ties introduced by the grid design. The final uncertainties result 
from the quadratic sum of these two uncertainty sources. 
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Fig. 3. Best fit model for the spectrum of HD 14818 degraded to R - 1000 and SNR - 100. The lines used in the analysis (marked 
in the plot) are well reproduced. The rest of transitions, some of them interstellar features, were not considered in the analysis. 



2.2.1. An example: HD 14818 

To illustrate the determination of the stellar parameters, we 
present the application of our algorithm to the case of the 
Galactic B-type supergiant HD 14818. A high quality optical 
spectrum is first degraded to the characteristics of our FORS2 
data, R = 1000 and SNR - 100. The;^'- calculations were ap- 
plied on individual wavelength windows, covering the spectral 
features shown in Table [T] Each wavelength range was carefully 
selected to include the entire line profile. Figure |2]displays the 
results obtained from the Monte Carlo simulation for each indi- 
vidual parameter As indicated, the spectral resolution, SNR and 
the set of lines were selected to represent the conditions of the 
spectra that will be analyzed in the forthcoming section (He is 
not included in our analysis since it is not available in the case 
of the NGC 55 stars). 

The quality of the final results is illustrated by Fig. |3] where 
the observed spectrum of HD 14818 is compared with a model 
atmosphere computed for the parameters obtained by our auto- 
matic analysis algorithm, desciibed in previous sections. Note 
that only the lines listed in Table [T] (marked in the figure) are 
considered and fitted for the analysis. 



2.2.2. Tests to the method 

Before applying the routine to the analysis of real data, we car- 
ried out a number of tests to check the reliability of the pro- 
posed methodology. In the first one we analyzed a sample of 
three synthetic spectra generated with the fastwind code. For 
the second test, we considered the spectra of three Galactic 
stars (HD 209975 09.5 lb, HD 38771 B0.5 la and HD 14818 
B2 la), originally from the lACOB spectroscopic database. In 
both cases, the test spectra were degraded to R - 1000 and 
SNR = 100 to simulate our FORS2 NGC 55 spectra. Taking 
into account that our FORS2 data do not include the Hff line, 
we decided to repeat the tests to evaluate the impact of includ- 
ing this line (columns 'Output' and 'Output+Ha' in Tables |2]-[3] 
respectively). 



Comparison to synthetic spectra 

Three sets of stellar parameters in the range of OB supergiant 
stars (log gf ~ 1 .0 - 1 .5 dex) were randomly chosen from our 
model grid and analyzed according to the method presented in 
Sect.|2] Table |2]presents the input parameters of the models and 
the stellar parameters recovered by the method. At this spectral 
resolution, different combinations of parameters, like microtur- 
bulence, silicon and helium abundance, can produce similar pro- 
files which has a clear impact on the uncertainties. Nevertheless, 
our algorithm recovers the input values within the errors. 

The analysis yields a better estimate of the wind parameter 
Q when Ha is incorporated. This also produces a slight variation 
in the rest of parameters, but always within the uncertainties of 
the previous results obtained without Ho-. It seems possible to 
constrain Q even without Ho-, just relying on the rest of Balmer 
lines (mainly H/3). Only for the very hot case (30900 K) we find 
large differences when Ha is not included in the analysis, with 
a shift of 0.20 dex respect to the input value. For the other two 
cases, the differences with respect to the input value do not ex- 
ceed 0.05 dex. 

Additional tests were performed with the spectra degraded to 
SNR - 50 (see Table |2]i. As could be expected, a wider set of 
models are compatible with the observations in this case, which 
translates in larger uncertainties. Nevertheless, the input param- 
eters and the recovered values are in good agreement, always 
within Icr. 



Comparison to Galactic stars 

The results obtained in the analysis of the three Galactic su- 
pergiants are collected in Table |3] al ong with t h e para meters de- 
rived by iRepolust et"ai1 (|2004 and lUrbaneial (|2004 . We have 
restricted the comparison to these two works because we want 
to minimize the possible effects introduced by the application 
of different stellar atmosphere codes; here, we are interested in 
the performance characteristics of our methodology. Both stud- 
ies used FASTWIND modcl s, although there are s ome unavoidable 
differences. For instance. IRepolust et al.l (|2004|) kept the 



Table 2. Results of the analysis of three synthetic spectra calculated with FASTWIND (parameters listed under column 'Input') degraded to R - 1000 and SNR - 100, 50. 





T309G329 


T265G280 


T134G166 




Input 


Output 


Output+Hff 


SNR = 5Q 


Input 


Output 


Output+Ho- 


SNR = 50 


Input 


Output 


Output+Hff 


SNR = 50 


reff(kK) 


30.9 


30.8 ± 0.8 


30.8 ± 1.1 


30.2 ± 2.0 


26.5 


26.1 ±0.8 


25.8 ± 1.1 


25.9 ± 1.6 


13.4 


13.3 ±0.6 


13.2 ±0.5 


13.0 ±0.7 


log g (dex) 


3.29 


3.25 ± 0.08 


3.26 ±0.08 


3.20 ±0.15 


2.80 


2.74 ±0.11 


2.72 ±0.11 


2.74 ±0.13 


1.66 


1.65 ±0.17 


1.66 ±0.07 


1.63 ±0.16 


4^(kms-') 


14.0 


13.6 ±3.3 


14.4 ±4.1 


13.0 ±6.4 


17.0 


17.5 ±4.0 


17.9 ±5.9 


15.0 ±8.0 


20.0 


22.5 ±5.3 


19.6 ±3.4 


20.5 ± 5.0 


He/H 


0.10 


0.12 ±0.03 


0.13 ±0.04 


0.12 ±0.05 


0.12 


0.15 ±0.04 


0.14 ±0.05 


0.13 ±0.05 


0.16 


0.15 ±0.04 


0.16 ±0.04 


0.12 ±0.05 


Si/Sie 


1.00 


1.01 ±0.25 


0.94 ± 0.24 


0.87 ± 0.28 


0.75 


0.78 ± 0.30 


0.66 ± 0.30 


0.82 ± 0.0.27 


0.80 


0.97 ± 0.23 


0.81 ±0.23 


0.89 ± 0.27 


loge 


-12.70 


-12.90 ±0.35 


-12.62 ±0.20 


-12.89 ±0.33 


-12.30 


-12.30 ±0.20 


-12.31 ±0.20 


-12.34 ±0.20 


-12.20 


-12.24 ±0.20 


-12.19 ±0.20 


-12.28 ±0.10 



Notes. The resulting parameters without Ha are gathered in column 'Output'. 
{'Output+Ha'). The SNR = 50 case considers Ha. 

Table 3. Results of the procedure for three Galactic stars degraded to R - 1000 and SNR - 100, 50. 



The process was also repeated including Ha, to estimate the effect of its absence on the stellar parameters 





HD 209975 


HD 38771 


HD 14818 




R04 


Output Output+Ho- 


SNR = 50 


U04 


Output 


Output+Ha 


SNR = 50 


U04 


Output Output+Ha 


SNR = 50 


reff(kK) 


32.0 


32.7 ±1.0 32.3 ±1.5 


32.0 ±3.3 


26.5 


27.2 ±1.6 


26.9 ± 2.0 


25.4 ±3.3 


20.1 


18.3 ±0.7 18.9 ±1.3 


18.9 ±1.9 


log g (dex) 


3.20 


3.36 ± 0.09 3.32 ± 0.08 


3.29 ±0.20 


2.90 


3.11 ±0.10 


3.11±0.11 


3.00 ± 0.20 


2.40 


2.51 ±0.11 2.45 ±0.14 


2.40 ±0.18 


4^(kms-l) 


10.0* 


14.7 ±5.3 14.4 ±5.0 


13.9 ±7.5 


18.0 


14.5 ±4.6 


14.1 ±4.3 


11.3 ±6.0 


17.0 


15.9 ±3.9 13.9 ±3.7 


11.9±5.9 


He/H 


0.10' 


0.11 ±0.04 0.11 ±0.04 


0.12 ±0.05 


0.10 


0.15 ±0.04 


0.16 ±0.04 


0.14 ±0.05 


0.15 


0.11 ±0.04 0.14 ±0.05 


0.12 ±0.05 


Si/Sie 


1.00* 


0.85 ±0.30 1.08 ±0.24 


0.83 ± 0.28 


1.00* 


0.91 ±0.29 


0.95 ± 0.26 


0.81 ±0.27 


1.00* 


0.87 ± 0.30 0.87 ± 0.29 


0.76 ± 0.27 


logG 


-12.67 


-12.69 ±0.38 -12.45 ±0.20 


-12.53 ±0.25 


-12.67 


-13.01 ±0.29 


-12.93 ±0.25 


-13.11 ±0.32 


-13.22 


-13.03 ±0.20 -13.15 ±0.20 


-13.18±0.11 



Notes . The results, summarized in this Table, with the errors ( 'Output') are compared with analyses performed (though with different observational data) bv lRepolust et al.l l l2004tl {R04) and lUrbaneial 
| |2004) {U04). The column 'Output+Ha' gives the values returned with the incorporation of Ha. The SNR = 50 case considers Ha. The parameters marked with '*' were kept fixed. 
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microturbulence fixed to 10 km s ', consi dering it as a seco ndary 
parameter; HD 14818 was analyzed by lUrbaneial (l2004l) with 
an early version of fastwind that did not include line block- 
ing/blanketing, which explains the difference in temperature 
(and gravity) with respect to our values. 

The main conclusion of these tests is that there is a good 
agreement between the values obtained by our algorithm applied 
to the low spectral resolution data and the previous studies based 
on high spectral resolution data. 

We note that there are important differences in log Q for the 
three stars, although the values are consistent within the errors. 
Synthetic Ha profiles depend not only on the wind parameters, 
but also on the effective temperature and surface gravity. Thus, 
changes in these two fundamental parameters will modify its 
shape with the consequent adjustment of the value recovered for 
Q. The uncertainties in the stellar parameters will propagate to 
log Q as well. The differences we have found could also be a re- 
flection of real changes in the observed profile, due to the use of 
Ha profiles collected in different observational campa igns (Ha 
is variable in these kind of stars, iMarkova et al.ll2005h . Finally, 
we should keep in mind that there could be (some) differences 
inherited from employing different versions of the same model 
atmosphere code. Nevertheless, in spite of these differences, the 
values recovered by our method are comparable with the ones 
suggested by these other works based on high spectral resolu- 
tion data. 

As in the case of the synthetic models discussed above, we 
performed the analysis of the observed spectra degraded also to 
SNR - 50. The outcome of this test is the same; the loss of 
information due to the lower SNR is reflected in larger uncer- 
tainties. 

Both sets of tests confirm the reliability of the technique 
in finding the main stellar parameters. The quality of the data 
(SNR) and the available observed wavelength range define the 
accuracy that can be achieved. We have shown that, by using 
enough information (see Table [TJ, it is possible to reach solid 
results, making our analysis algorithm a very promising tool for 
the analysis of large collections of optical spectra of OB stars. 

Our main focus is the analysis of supergiant stars. Hence we 
have designed the grid and selected the spectral features for the 
analysis accordingly. However, this methodology can be applied 
to other spectral types/luminosity classes. Of course, this would 
require different diagnostic transitions, a different set of models, 
and would present d ifferent challenges. The reader is referred to 
iLefever et al.l (1201 Oh for a thorough discussion and application 
of a very similar methodology in the case of Galactic dwarf and 
giant stars. 

2.3. Abundance determination 

Once the fundamental parameters have been determined, we pro- 
ce ed with the ana l ysis of the ch emical abundances. As shown 
by lUrbaneiaet al.l (l2003l l2005ah . it is possible to derive indi- 
vidual chemical abundances from low resolution optical spec- 
tra of early B-type supergiants because of the relative low num- 
ber density of metal lines (hence minimizing the blends of lines 
from different species) and because there is a good number of 
relatively strong features that can be easily detected. Following 
these considerations, the methodology applied in the chemical 
abundance analysis relies on simultaneously modeling all the di- 
agnostic lines for a given species. 

Using the derived stellar parameters, a new set of tailored 
models is computed for each star under analysis by varying the 
abundances of the relevant species (C , N , O and Si ) in steps 



of 0.20 dex. The chemical analysis is performed in two comple- 
mentary steps. First, an automatic x^ fitting algorithm is used. 
However, the main features were sometimes weak or blended, 
and a second visual check is required. 

The abundance uncertainties are estimated from this new set 
of models, taking into consideration the SNR of the observed 
spectra. Therefore, the range of abundances defined by the un- 
certainties, for each individual species, account for the feature- 
to-feature scattering, in a similar way as it would be done in a 
classic analysis. 

The precision achieved in the chemical analysis depends not 
only on the quality of the spectra, but also on the spectral type 
and on the reliability of the atomic models. For a mid B-type, 
most of the considered transitions cannot be detected (either they 
are not present at these temperatures, or they are too weak to be 
detected at this low resolution and SNR). Note, however, that the 
Mg II and Si II features are stronger for late B-type stars, so we 
can extract accurate information for those species. In a similar 
fashion, for a given spectral type and SNR, the abundance un- 
certainties will depend on the metallicity, with the expectation 
that errors become larger with decreasing metallicity, eventually 
reaching the limit where only upper limits could be placed. 

For each individual element, we attended to a particu- 
lar group of spec tral features (see ;Urbaneia et al. 2005al and 
lEvans et aDl2007h . Briefly, the lines considered for the chemi- 
cal analysis, and summarized in Table |4] were: 

- Silicon. We used the transitions of Si III 4552 - 67 - 74 A, 
Si II 4128 -30 A and Si IV 41 16 A. The lines of Si II are only 
available at low temperatures, B2 and later spectral types. On 
the other hand, the profile of Si IV will be accessible at late 
O and early B-type stars. 

- Oxygen. The analysis looks for the best fit to the transitions 
around 4076, 4319, 4350 and 4416 A; at low spectral resolu- 
tion they are in fact blends. There are additional O II features 
around 4585 and 4700 A, but they are blended with other el- 
ements. The latter will be used as a consistency check for the 
abundance derived using the other O II lines. 

- Nitrogen. Some of the most prominent N features in the 
wavelength range covered by our FORS2 spectra are blended 
with other elements (for example around 4650 A), but there 
are still some isolated features that can be used to constrain 
the nitrogen abundance. The study was centered on 3995, 
5007 and 5045 A. Note that if nebular subtraction is not ac- 
curate, the [O III] 5007 A transition could severely affect the 
overlapping nitrogen line. 

- Carbon. Given the spectral quality and the wavelength 
range (~ 3900 - 5000 A) the strongest CII transition is 
4267 A. We consider that the results based (only) on this 
line are currently not as reliable as for the rest of the species 
(see lLennon et al.ll2003l) . Alternative lines are 3919 - 21 A, 
though too weak in many cases. The C III transitions around 
4650 A, while blended with O and N, could serve as a sec- 
ondary check. 

- Magnesium. Its abundance is determined using the line of 
Mg II 448 1 A. Note that this transition could be blended with 
Al III 4479 A dLvubimkov et alJl2005h . 

Figure |4] illustrates the results of the chemical analysis of 
HD 14818. The main diagnostic lines used in the elemental 
abundance determination, identified in the figure, show a good 
match with the final model. 
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Fig. 4. Comparison of the final fastwind model (red), including the different species considered in the chemical analysis (He, C, N, 
O, Mg and Si) with the low resolution spectrum of HD 14818 (R - 1000 and SNR - 100) in black. The plot also displays the best 
model -0.20 and +0.20 dex in the abundances of the elements included in blue and green respectively. The final model provides a 
good match to the observation. The key diagnostic metal features are identified. 
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Table 4. Metal features considered for the chemical analysis. 



Silicon 


Carbon 


Nitrogen 


Oxygen 


Magnesium 


Ion 


^(A) 


Ion 


A (A) 


Ion A (A) 


Ion 


^(A) 


Ion A (A) 


Si IV 


4116 


CII 


3919 


Nil 3995 


on 


4076 


Mgll 4481 


Sim 


4552 


CII 


3921 


Nil 5007 


on 


4319 




Si III 


4567 


CII 


4267 


Nil 5045 


on 


4350 




Si III 


4574 








on 


4416 




Sill 


4128 














Sill 


4130 
















Fig. 5. Distribution of the objects selected for the analysis 
over NGC 55. The image was taken from the DSS archive 
( http://archive.stsci.edu /cgi-bin/dss Jbrm^ . The black cross marks 
the optical galactic center. 



3. Quantitative spectral analysis of NGC 55 B-type 
supergiant stars 

In COS we presented the first spectral catalog of massive 
blue stars in NGC 55. Very briefly, optical spectra of ~200 
sources were collected wi th the FOcal Reducer/lo w dispersion 
Spectrograph 2 (FORS2, lAppenzeller et al.lll998l) at the Very 
Large Telescope (VLT-UT2). The instrument was equipped with 
the 600B grism, providing optical spectra in the range ~ 3900 - 
6000 A at a resolving power R ~ 1000. A complete description 
of the data and the reduction process can be found in the afore- 
mentioned reference. 

From this previous work, we selected twelve B-type super- 
giants for a detailed quantitative analysis. The selection was 
based on the good SNR, spectral type and the lack of obvious 
contamination by other sources, including strong nebular lines 
(when possible). The selection of spectral types between late O- 
and early B-types guarantee that the main spectral features re- 
quired for the analysis (see Table [T]! are available. The stars are 
spatially distributed across the galaxy (see Fig. |5]l, which will 
allow us to investigate the distribution of its chemical composi- 
tion. Table |5] summarizes all the relevant information, as well as 
provides revised photometry. The stars are identified following 
COS. 



3.1. Stellar parameters 

The stellar parameters of our sample of 12 NGC 55 B-type su- 
pergiants, obtained from the application of the analysis method- 
ology described in the previous sections, are listed in Table |6] 
The left side of Figs. |6] and [T] show the comparison of the ob- 
servations with tailored models computed for the parameters de- 
rived in this work. It can be seen that these final models provide 
a good match to the (in some cases rather noisy) observed spec- 
tra. The right side of the figures display log x' isocontours on 
the Teff - log g plane. Each black dot represents a model in the 
grid, whilst the white dots identify those models fulfilling the 
y~ - y^ ■ criteria defined in Sect. 12.21 

^ ^ min ' ' 

The impact that the SNR has on the derived parameters (un- 
certainties) can be gauged from these plots. For stars like A_l 1 
and B_31, the models that reproduce the observations enclose a 
relatively smaller area in the Tefj - log g plane than in the case 
of Cl_45. The higher quality is clearly reflected in the range of 
stellar parameters, i.e. synthetic models, that are compatible with 
the observed spectra. Note how the wide range of temperatures 
and gravities covered by our grid prevented, for our 12 stars, the 
presence of border-effects, i.e. objects located too close to the 
limits of the model grid. 



3.2. Chemical abundances 

With the stellar parameters in hand, we performed the analysis of 
the elemental chemical abundances as described in Sect. 12. 31 The 
results are compiled in Table |7] The characteristic metallicity, 
given in the last column, is obtained averaging the differences 
of the derived O, M g and Si abundances relative to the solar 
reference, taken from lAsplund et aD(l2009i) . Projected distances 
to the galactic center are given in the second column, in units of 
the semi-major axis. 

The final models presented in Figures |6]and|7]include the de- 
rived elemental abundances. Overall, these figures show a good 
agreement between the final tailored models, computed for the 
final abundances, and the observed spectra. The average abun- 
dance uncertainty in these chemical analyses is ~0.25 dex, with 
the case of Cl_45, with its poor SNR, reaching 0.30 dex. Our re- 
sults indicate that it is possible to chemically characterize these 
stars with good precision, even for cases of relatively low SNR, 
suchasCl.53orCl.45. 

We note here that the derived Mg abundances, based solely 
on the Mgll 44S1 A feature, could be affected by a blend 
with Aim 4479 A. Fol lowi ng the arguments in t he works by 
iLvubimkov et all (|2005|) and iDunstall et all ( 1201 lb . a small ef- 
fect woul d be expected, given the low metallicity of NGC 55. 
Moreover. lUrbanei a et al . (2005b) suggested that there is a neg- 
ative luminosity dependence of the strength of the Al III 4479 A 
line, making its contribution even less relevant for supergiant 



10 



N. Castro et al.: Quantitative Study of Massive Blue Stars in NGC 55 




4000 4200 4400 4600 

Wavelength (A) 



4800 



5000 



10 15 20 25 30 
T^f (kK) 



Fig. 6. Left panels: NGC 55 spectra (black) and best fit model (red) obtained by our automatic method plus the detailed abundance 
analysis. The lines used in the stellar parameter and chemical abundance determination are marked in the top-left chart. Right panel: 
log x^ difference between observations and the entire grid models in the Teff -log g parameter space. The^^ difference is quantified 
with the color, according to the scale shown at the right of each plot. The black dashed lines mark the model with minimum ;if^. 
Black dots represent all the grid models; the models whose parameters have been averaged to produce the stellar parameters are 
highlighted with white dots. 
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Table 5. NGC 55 stars analyzed in this work. 



ID 


RA(J2000) 


DEC(J2000) 


SpT 


SNR 


V 


I 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


A.8 


0:15:44.18 


-39:14:58.20 


09.71 


62 


20.067 ± 0.005 


20.079 ± 0.010 


C1.44 


0:15:11.29 


-39:12:34.55 


BOI 


79 


19.360 ± 0.005 


19.331 ±0.007 


C1.9 


0:14:45.54 


-39:12:38.17 


BlI 


86 


19.304 ±0.004 


19.504 ±0.009 


C1.13 


0:14:51.86 


-39:10:57.35 


BlI 


60 


19.222 ± 0.004 


19.254 ±0.006 


C1.45 


0:15:11.31 


-39:12:50.40 


BlI 


38 


19.744 ± 0.007 


19.714 ±0.012 


A.17 


0:15:51.42 


-39:15:57.60 


BlI 


79 


20.426 ± 0.007 


20.157 ±0.008 


D.27 


0:14:28.73 


-39:10:17.76 


B2I 


76 


19.395 ±0.004 


19.373 ±0.007 


A.27 


0:15:58.03 


-39:14:27.60 


B2I 


66 


20.140 ±0.006 


20.050 ± 0.010 


C1.53 


0:15:14.64 


-39:13:36.12 


B2.5I 


45 


20.174 ±0.005 


19.837 ±0.007 


B.31 


0:15:32.22 


-39:14:40.20 


B2.5I 


95 


19.555 ± 0.005 


19.561 ±0.007 


A.26 


0:15:55.45 


-39:15:30.60 


B2.5I 


96 


19.584 ±0.005 


19.197 ±0.005 


A.11 


0:15:45.17 


-39:15:56.15 


B5I 


130 


19.005 ± 0.004 


18.874 ±0.005 



Notes. The columns list: (1) star identification; (2) right ascension (hh:mm:ss); (3) declination (dd:mm:ss); (4) spectral type; (5) signal-to-noise ra- 
tio of the spectrum; (6) apparent V magnitude; (7) apparent I magnitude. New calibrations have reveled issues in the published NGC 55 photometry 
catalog in COS, here a revised photometry is shown for the stars targeted (Castro et al. in prep). 

Table 6. Stellar parameters obtained for the NGC 55 stars by the automatic algorithm presented in this work. 



ID 


Teff 


logs 


^ 


He/H 


logg 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


A.8 


27700 ± 1200 


2.91 ±0.13 


20.3 ±5.8 


0.14 ±0.05 


-12.79 ±0.33 


C1.44 


26200 ± 1500 


2.83 ±0.15 


17.8 ±6.3 


0.10 ±0.05 


-12.93 ±0.30 


C1.9 


17400 ± 1500 


2.37 ±0.18 


11.8±3.3 


0.15 ±0.05 


-13.37 ±0.27 


C1.13 


24200 ± 1400 


2.69 ±0.15 


17.7 ±5.7 


0.15 ±0.05 


-12.60 ±0.20 


C1.45 


21500 ±2500 


2.73 ± 0.24 


19.0 ±6.1 


0.12 ±0.05 


-13.65 ±0.36 


A_17 


22500 ± 2300 


2.85 ±0.16 


20.7 ±5.3 


0.10 ±0.04 


-13.90 ±0.33 


D_27 


19400 ± 1600 


2.41 ±0.15 


14.4 ±4.6 


0.10 ±0.05 


-13.18 ±0.20 


A_27 


16300 ± 1100 


2.13 ±0.19 


20.4 ± 5.7 


0.12 ±0.05 


-12.90 ±0.20 


C1.53 


18200 ±2000 


2.54 ± 0.26 


16.4 ±5.8 


0.10 ±0.05 


-13.73 ±0.35 


B.31 


16700 ± 900 


2.09 ±0.13 


20.4 ± 5.6 


0.14 ±0.05 


-12.83 ±0.20 


A_26 


16100 ± 1000 


2.15 ±0.17 


17.1 ±6.5 


0.12 ±0.05 


-12.86 ±0.20 


A_ll 


14400 ± 1000 


1.94 ±0.18 


11.8 ±3.7 


0.10 ±0.05 


-12.64 ±0.20 



Notes. The columns list: (I) star identification; (2) effective temperature (K); (3) surface gravity (dex); (4) microturbulence (kms 1); (5) helium 
abundance (by number); and (6) the Q parameter. 



stars, even more since Mg II 448 1 A presents the opposite behav- 
ior, with the feature strengthening with increasing luminosity. 

3.3. Comments on individual targets 

Here we discuss some details of individual targets, along with 
problems encountered during their analysis. 

- A_8. This is an 09.7 I star whose main features are well 
represented by our final model. At this hot temperature 
(27700 K) the Mg II line has vanished and the Si III lines 
are weak. At the same time. He II lines are present, allow- 
ing for a precise determination of the effective temperature. 
The Si IV line, blended with O II, is well reproduced. The 
apparent slight discrepancy around the O II 4079 A line is 
most likely due to an effect of the normalization inside the 
U6 wing. Note also that the O II 4414 - 16 A blend is weaker 
than the prediction. Hell 4686 A shows also a clear mis- 
match but, without a reliable wind estimation we cannot go 
further in its analysis. 

- Cl_44. The main spectral features are well reproduced by the 
final model, with the exception of some transitions like: He II 
4686 A, possibly affected by the wind, and not considered 
in the analysis; or Si IV 4089 A, not included because of its 



blend with O II. The O II doublet at 4414- 16 A is also poorly 
reproduced, though the rest of the oxygen lines are consistent 
with the obtained abundance. At this temperature the Mg II 
line is quite weak. The transition of N II 5007 A is filled by 
nebular contamination. 

Cl_9. The final model shows a good global match to the ma- 
jority of the considered transitions. The lines of O, Mg and 
Si are not particularly strong at this temperature, but they are 
well reproduced. The temperature derived corresponds to a 
spectral type cooler than the one assigned. 
Cl_13. This Bl supergiant presents a very nice match be- 
tween the final model and the observed spectrum. The blends 
of on 4079 A and 4319 A show the largest difference, but 
still within the abundance uncertainties of +0.25 dex. Some 
contamination by nebular emission is still apparent in the 
N II 5007 A line. The Mg II is weak, as expected for this 
temperature. 

Cl_45. This spectrum has the lowest SNR ratio in the sam- 
ple and this is clearly reflected in the errors derived for the 
stellar parameters. Figure |6]shows a large spread of compat- 
ible models in the Tgff - log g plane. Nonetheless, the model 
corresponding to our solution reproduces well the main fea- 
tures. The observed N II 3995 A line is not well represented 
by the model, and the transition at 5007 A is contaminated by 
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Fig. 7. Figure |6] continued. 



nebular emission. As typical for a Bl I, the Mgll transition 
is weak. 



- A_17. This object shows a discrepancy at the continuum 
around H6, likely due to the rectification around H6 wings. 
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Table 7. 

O. 



Stellar abundances determined for Si , Mg , C , N and 



ID 


PlPo 


esi" 


ec" 


en" 


fo" 


fMa" 


[Z/Zo] 


O 




7.51 


8.43 


7.83 


8.69 


7.60 




A.8 


0.64 


7.18 


7.27 


7.33 


8.28 


— 


-0.37 


CI. 44 


0.22 


7.08 


7.71 


7.56 


8.07 


— 


-0.53 


CU 


-0.11 


6.80 


7.23 


7.95 


8.66 


6.80 


-0.51 


C1.13 


-0.06 


7.07 


7.60 


7.63 


8.45 


6.92 


-0.45 


C1.45 


0.22 


7.13 


8.07 


8.17 


8.34 


7.00 


-0.44 


A.17 


0.75 


7.06 


7.77 


7.64 


8.37 


7.02 


-0.45 


D.27 


-0.32 


7.01 


7.37 


8.28 


8.26 


7.27 


-0.42 


A.27 


0.80 


7.04 


7.77 


8.52 


8.07 


6.98 


-0.57 


C1.53 


0.28 


7.25 


7.52 


8.22 


8.63 


7.28 


-0.21 


B.31 


0.50 


7.10 


7.22 


8.25 


8.59 


7.26 


-0.28 


A.26 


0.78 


7.43 


7.48 


8.22 


8.60 


7.33 


-0.15 


A.11 


0.68 


7.17 


— 


— 


— 


7.20 


-0.37 
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Notes. The projected radial distance of the stars to the optical center of 
NGC 55 normalized by the major semi-axis (p„ ~ 16') is displayed as 
well (plpo)- Negative values of pi po mean that these stars are located 
westward of the galactic c enter. The solar abun dances used as a refer- 
ences (O) were taken from lAsplund et aLl(l2009l) . 



(ii) 



ex = log (X /H ) -I- 12 (by number) 



Fig. 8. The sample stars in the HR diagram. D a shed L ines mark 
the evolutionary tracks of iMaeder & Mevnell (1200 lb at SMC 
metallicity with an initial rotational velocity of 300 km s"'. 
Those stars that have shown large discrepancies between the 
spectroscopic mass and the evolutionary one are marked by gray 
squares or diamonds (see text for additional details). 



The parameters and abundances derived reproduce well the 
rest of the spectrum. The magnesium transition is weak. Note 
in Fig.|6]the lack of symmetry in;^'- distribution around the 
average values. There is a plume of suitable models towards 
high temperatures that also fulfill our goodness-of-fit crite- 
ria. 

- D.27. The line of He I at 4144 A presents a mismatch. Note 
that this line is not included in the analysis. 

- A.27. This B2 I is well represented by the final model. 
Besides from the mismatch around the 5007 A region, the 
rest of the spectral features are well reproduced. 

- C1.53. In spite of the low SNR spectrum obtained for this 
star, the match between the main features analyzed and the 
best model is quite good. The spectrum also displays some 
residuals of nebular contamination at [O III] 5007 A. 

- B.31. This B2.5 supergiant star is nicely reproduced by our 
model, as it is shown in Fig.|7] The main transitions used in 
the analysis are well reproduced. 

- A.26. This B2.5 I star shows a good agreement with the 
model. Its better SNR ratio results in a well constrained set of 
parameters. The oxygen lines are weaker at this temperature 
but they are all well modeled except for O II 4079 A. The 
latter issue may be caused by the continuum rectification, as 
in A.17, or residuals from cosmic ray subtraction. 

- A.ll. This star has the best SNR spectrum, and this is re- 
flected in the small uncertainties. Its effective temperature 
suggests a spectral type a bit earlier than the one proposed 
by COS. At these cool temperatures the transitions of Oil, 
N II and C II are very weak or have even completely van- 
ished. The Si III lines are weak but the transitions of Si II 
(although blended with O II lines) become stronger and al- 
low us to derive effective temperature and constrain the sili- 
con abundance. 

4. Discussion 

In this section we discuss the results obtained for our sample of 
12 B-type supergiant stars in the galaxy NGC 55. The physical 
characterization provided by the stellar parameters and chemi- 



cal abundances, supplies us with the necessary information to 
discuss their evolutionary status, as well as to carry out a com- 
parison with the predictions of current evolutionary models. We 
adop t a distance rnodulus to NGC 55 of /i = 26.434 + 0.037 mag 
from lGieren et al.l (l2008l) . 



4. 1 . Stellar properties 

Table |8] gathers the fundamental stellar properties derived for 
the stars. The color excess of each object was calculated us- 
ing observed photometry (see Table |5) and the synthetic colors 
obtained from the final tailore d models, adopting the extinction 
curve bv lCardelli et al.l (Il989h and a total-to-selective-extinction 
ratio Rv = 3.1, although several authors have sho wn that high Ry 
are n ot rare for massive stars (se e, for instance, iBonanos et al.l 
I2OIII or lBestenlehner et aDl20I ih . Only for one of the objects, 
C1.9, we find a non-physical value (i.e. negative), although very 
small and compatible with zero within the error bars. The high 
inclination of the galaxy, as well as the fact that these objects 
could be, to some extent, surrounded by ionized gas, could 
bias the observed photometry which would explain thi s nega- 
tive value, an effect also suggested bv lEvans et al.l ( l2007h . Three 
other objects, A.17, A.26 and C.53, present high reddening val- 
ues, ~0.4 mag, whilst the rest of the sample show E{B-V) value s 
consistent with the mean value derived bv iGieren et al.l (|2008|) , 
E(B -V) = 0.127 + 0.019 mag, from multi-wavelength observa- 
tio ns of Cepheid stars . 

IWebster & Smiflil (Il983h studied a sample of 7 H II regions 
distributed along the southern half of NGC 55, approximately 
covering the same range in galactocentric distances as our stel- 
lar sample. From their published C {Hfi) values we can obtain 
redde ning values by adopting E{B - V) = 0.676 C(H/3) dSeatonl 
II979I) . The HII regions present reddening in the range 0.16- 
0.32 mag, with a simple mean of E(B - V) = 0.24 + 0.07 mag, 
with 3 H II regions showing E{B - V) > 0.3 mag. The mean 
value of our sample, E{B -V) = 0.20 + 0.14 mag, compares well 
with the nebular mean, albeit with a larger scatter. 

The distribution of the stars in the HR diagram is 
shown in Fig. [8] together with the evolutionary tracks by 
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Table 8. Photometry and fundamental parameters recovered from the sample analyzed in NGC 55. 



ID 


logg/^ 


Mbol 


E{B - V) 


BC 


log L/Le 


R/Re 


Mspcc-/A^o 


Mevoi/Mo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


A.8 


1.14 ±0.07 


-9.60 ±0.11 


0.177 ±0.014 


-2.686 ± 0.096 


5.74 ± 0.04 


32.2 ± 1.1 


30.7 ± 10.3 


32.7 ± 3.5 


CI. 44 


1.16 ±0.07 


-10.24 ±0.11 


0.199 ±0.017 


-2.547 ± 0.087 


5.99 ± 0.04 


48.2 ±3.1 


57.4 ±23.5 


46.8 ±5.8 


C1.9 


141 ±0.07 


-8.57 ±0.13 


-0.049 ± 0.016 


-1.592 ±0.119 


5.32 ± 0.05 


50.7 ± 5.6 


22.0 ±11.6 


21.7 ±1.8 


C1.13 


1.15 ±0.09 


-9.98 ±0.15 


0.126 ±0.020 


-2.375 ±0.130 


5.89 ± 0.06 


50.1 ±2.4 


44.9 ± 17.6 


40.1 ±4.7 


CI. 45 


1.40 ±0.12 


-9.26 ± 0.24 


0.171 ±0.024 


-2.041 ± 0.222 


5.60 ± 0.09 


45.6 ±5.7 


40.9 ± 27.7 


28.5 ±4.0 


A.17 


1.44 ±0.08 


-9.36 ± 0.26 


0.377 ± 0.026 


-2.185 ±0.249 


5.64 ±0.11 


43.6 ±3.7 


49.3 ± 22.3 


29.8 ± 4.3 


D.27 


1.26 ±0.10 


-9.29 ± 0.22 


0.137 ±0.021 


-1.824 ±0.206 


5.61 ± 0.09 


56.7 ± 3.7 


30.2 ± 12.4 


28.7 ±3.5 


A.27 


1.28 ±0.09 


-8.16 ±0.13 


0.151 ±0.019 


-1.393 ±0.114 


5.16 ±0.05 


47.7 ± 3.5 


11.2 ±5.7 


18.7 ± 1.6 


C1.53 


1.50±0.11 


-9.13 ±0.21 


0.394 ± 0.021 


-1.653 ±0.195 


5.55 ± 0.08 


60.1 ±7.5 


45.7 ±33.1 


26.0 ± 3.0 


B.31 


1.20 ±0.08 


-8.58 ±0.12 


0.078 ± 0.018 


-1.460 ±0.096 


5.33 ± 0.05 


55.3 ±3.0 


13.7 ±4.9 


21.1 ± 1.6 


A.26 


1.32 ±0.09 


-9.44 ±0.13 


0.392 ± 0.017 


-1.370±0.111 


5.67 ± 0.05 


88.2 ±5.8 


40.1 ± 18.3 


29.2 ± 3.0 


A.11 


1.31 ±0.08 


-9.07 ±0.12 


0.153 ±0.014 


-1.164±0.111 


5.52 ± 0.05 


93.0 ±7.6 


27.5 ± 13.7 


25.9 ±2.2 



Notes. The columns list: (1) star identification; (2) flux-weighted gravity (dex); (3) bolometric magnitude; (4) color excess in (B-V); (5) bolometric 
correction; (6) stellar luminosity; (7) radius; (8) spectroscopic mass; and (9) evolutionary mass. 
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Fig. 9. Left: Relationship between the spectroscopic and evolutionary masses derived in NGC 55 (the one-to-one trend is marked 
with a black solid line). Right: The logarithm ra tio of masses versus th e ir luminosity in this work (black dots), to gether with the 
results published previously in other galaxies by iKudritzki et alj ( l2008h . lUrbaneia et al.l ( 120081) and lU et all ( |2009|) (open circles). 
Those objects that show higher evolutionary mass compared with the spectroscopic counterpart are highUghted with gray squares 
for quick contrast with Fig. |8] The opposite behavior is indicated with gray diamonds. 



iMaeder & Mevnell (1200 ih . According to the metallicity of the 
sample (Table |7]i, we consider a linear interpolation between 
evolutionar y tracks com puted forjhe Small Magellanic Cloud 
metalhcitv dMaeder & M evnet 200 1!) and solar metallicity tracks 
dMevnet & Maeden I2003D . The evolutionary masses derived 
from the interpolation of these tracks are shown in the last col- 
umn of Table [8] These results were also checked wit h the recent 
evolut ionary calculations for LMC metallicity by iBrott et al] 
(1201 ih . finding very similar results. 

The comparison of spectroscopic and evolutionary masses 
has been an important source of discrepancy between the evolu- 
tionary theories and t he stellar atmosphere modeling for decades 
(iHerrero et al.l Il992h . Improvements in both fields have mini- 
mized this issue. Nonetheless, a systematic shift can still be mea- 
sured in the analysis of B -type supergian t stars. For instance, 
the analysis carried out bv lU et akl (l2009h in M 33 revealed an 
average difference of 0.06 dex between spectroscopic and evo- 
lutionary masses. Although small, it is still a systematic issue. 
Figure |9] displays the relationship between these two mass es- 
timations for our NGC 55 B-type supergiant stars. The errors 
in the spectroscopic masses are large enough to make the full 



sample consistent with the one-to-one relation (left panel in Fig. 
|9]i; half of the sample shows a very good agreement between 
both measurements and no systematic trend is evident, whilst 
four stars (A.17, A.26, CI. 45 and CI. 53, marked in the figure 
with gray squares) show significantly higher evolutionary than 
spectroscopic masses. Table [8] reveals that three of them show 
the highest color excess in the sample, which could point to a 
nebular contamination effect in the observed photometry. On the 
other hand, A.27 and B.31 show the opposite behavior (marked 
with gray diamonds in Fig. |9]l. Due to the spatial and spectro- 
scopic resolution we cannot rule out additional unresolved com- 
panion(s) that would be affecting the photometry (but are not 
evident in the optical spectra). 



4.2. Flux-weighted Gravity-Luminosity Relationsliip 

Figure [To] shows the flux weighted gravity-luminosity relation- 
shi p for the NGC 5 5 star s, together with the resu lts published 
by IKudritzki et alj JIOOSl) . lUrbaneia et all dloosl) and lU et all 
(l2009h . The location of the NGC 55 stars shows a good agree- 
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Table 9. Element abundances of Si, 
averaged over the whole sample. 



O , Mg and estimated Z 



2.4 2.2 2.0 1.8 1.6 1.4 1.2 1.0 
log gp (dex) 

Fig. 10. Position for our NGC 55 sam ple in the FGLR (blac k 
dots). Previous studi es pe rformed by Kudritzki et all (l2008h . 
lUrbaneia et"ai1 (|2008|) and lU et al.1 (|2009|) in different galaxies 
ar e also displayed (open c ircles), the linear relationship proposed 
bv lKudritzki et alJ (l2008h (see Eq.[T]i is drawn with a solid black 
line. The NGC 55 stars that have shown stronger discrepancies 
between the spectroscopic and evolutionary masses are marked 
as in Fig. [8](see Sect. l4.1l for additional details). 



ment with these studies, worst for those stars for which we 
have found disagreements between the different mass estima- 
tions. The other six stars follow the same trend with a slight 
shift towards higher bolometric magnitudes, alt hough they are 
within the observed scatter of the distribution (iKudritzki et al.l 
|2008|) . An independent determination of the distance to NGC 55 
based on the FGLR is deferred to a future publication, since in- 
formation from BA supergiant stars, with log gf > 1.5 dex, is 
mandatory to properly determine the distance modulus. 

4.3. Evolutionary chemical status 

Current massive star evolutionary models, accounting for the ef- 
fects of mass-loss and rotation, predict a tight relationship be- 
tween N/C and N/O surface ratios as a consequence of the 
mixing with CNO processed material from inner layers during 
the stellar evolution (e.g. iChiosi & Maeder 1986) . Previous stud- 
ies have observati onally found this rela tionship in Galactic stars 
(see for instance iPrzvbilla et al.l 1201 Ol and references therein), 
although the large uncertainties leave open a broad range of 
interpr etations. The detailed study presented by IPrzvbilla et al.l 
( 1201 Ol) on a sample of Galactic B-type dwarfs and A-type su- 
pergiant stars of the solar neighborhood revealed a very good 
agreement with the theoretical predictions of single star evolu- 
tionary models, in the particular range of stellar masses sampled 
by these objects, -20-40 M©. It has to be still proven that this 
is also the case for other mass ranges. Our sample in NGC 55 
contains an heterogeneous group of stars throughout the disc of 
NGC 55, nonetheless the left panel of Fig. (TTj shows that these 
stars follow the theoretical predictions in a qualitative sense. 

It would seem that our derived N abundances define two 
groups of stars, with the first one clustering around 7.62+0.22 
dex (simple mean and standard deviation) and a second one at 
8. 28 ±0.12 5 dex. A dopting as N baseline the HII region values 
derived bv lWebster & Sm ith (1 983i ), 6.63+0. 10 dex and a mean 
N/O - 0.015, all our stars show a high degree of N processing: 
N/O - 0.15 and 0.69 (mean values for each group). Interestingly 



[X/H] 


Mean ± cr 


a 


b 


Si 


-0.40 + 0.15 


0.05 + 0.07 


-0.42 + 0.03 





-0.30 + 0.21 


-0.06 + 0.18 


-0.28 + 0.09 


Mg 


-0.49 + 0.18 


0.11+0.15 


-0.53 + 0.08 


Z 


-0.40 + 0.13 


0.09 + 0.10 


-0.43 + 0.05 



Notes. The columns a and b provide the coefficients of the linear regres- 
sion (see Fig. ll2t . which accounts for the spatial chemical distribution. 



enough, similar values of enrichment have been recently re- 
ported for O type stars in the LMC by iRivero Gonzalez et al.l 
(I2OI2I) . These authors find two distinct groups of strongly en- 
riched objects, with N abundances 7.5 dex and 8.1 dex, with the 
LMC baseline N abundance being at 6.9 dex. Our B-type su- 
pergiant stars in NGC 55 show a remarkable agreement with the 
N abundances of O type stars in the LMC. This strongly sup- 
ports the idea that our B supergiant stars belong to a young pop- 
ulation evolving away from the Main Sequence towards the red 
part of the HR diagram. We see no indication of any object be- 
ing in a blue-loop, i.e. being a post Red Supergiant object: be- 
sides the previous discussion, current evolutionary models pre- 
dict that the blue loops cannot reach the temperatures of our ob- 
jects. Therefore, the main conclusion is then that none of these 
objects is in an advanced evolutionary stage. 

4.4. Metallicity distribution in the disk of NGC 55 

Previous works based on the study of the emission line spectra of 
H II regions have found that the present-day m etallicity (oxygen) 
of NGC 55 is very si milar to that of the LMC (IWebster & Smithl 
ll983l:lDavidgdl2005b . From our sample of 12 B-type supergiant 
stars we find a mean metallicity of -0.40 ±0.13 dex (see the first 
column in Table |9li, a value quite close to the LMC metallicity 
( jHunter et al. 2007 ). We reached a similar conclusion with our 
previous qualitative analysis (C08). 

We also calcul ated the O/H abundance ratio for the four H II 
regions for which IWebster & Smithl (Il983h reported the detec- 
tion of the [O III] 4363 A auroral line. This allows us to provide 
'direct' nebular abundances, which are independent of various 
calibrations of strong-line methods. Electron temperatures were 
calculated with the temden program in IRAFQ using the [O III] 
4363 /(4959 + 5007) line r atio and updated atomic parameters, 
as in lBresolin et al.l ( |2009|) . The O^ and O^^ ionic abundances 
were then calculated with the program ionic. We then obtained 
O/H as result of the sum of these two ionic abundances. Figure 
ITS] displays the excellent agreement between these new abun- 
dances (shown as stars in the figure) and the ones derived from 
our sample of B-type supergiants. 

The inclination of NGC 55, along with its apparently irreg- 
ular shape, make this galaxy a difficult object for morphologi- 
cal classification. The presence of a metallicity gradient across 
the galaxy would hint to a possible spiral disc. The study of 
IWebster & Smithl (Il983h did not find any trace of spatial varia- 
tions in the southern part of the galaxy. To investigate this issue, 
the distribution of silicon, oxygen and magnesium were mea- 
sured from our sample of B-type supergiants. We show the spa- 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 11. Left panel: [N/C] versus [N/O] ratios (by number) relative to solar values derived in our NGC 55 sample, the average 
uncertainty is plotted in th e right-bottom comer The evolutionary tracks for three different m asses (20, 40 and 60 Mp), are displayed 
for solar (solid black lines. iMevnet & Maedeil2003h and SMC metallicity (dashed gray lines. lMaeder & Mevnetl200lh . Right panel: 
position of the sample stars in the [N/O ] - Teff diagram, using the same color and symbol code that in the left panel. 



tial trends in Fig. [12] The individual elemental abundances hint 
to an almost null gradient. In order to quantify these results, 
we fit radially dependent gradients, [X/H ] = a {plpo) + b, to 
our stellar data. We exclude from the fits those objects with val- 
ues beyond ±2 cr of the mean value. In the previous expression, 
[X/H] = log (X/H) - log (X/H )o, XI H represents the abun- 
dance of each element relative to H by number, {X/H)q is the so- 
lar reference, and p is the projected galactocentric distance (the 
semi-major axis, p,, ~ 16' for NGC 55). The parameters of the 
linear regression, a and b are collected in Table |9] Considering 
the errors in the regression coefficients (Table|9]l, all the elements 
show a spatial di stribution consistent with no gradient, support- 
ing the results bv lWebster & Smithl(ll983h . but in this case based 
not only O, but also Mg and Si. Nonetheless, it would be highly 
desirable to expand the sample in the western half of the galaxy, 
before drawing a definitive conclusion. The right side of Fig. [12] 
shows the 2D abundance distribution over NGC 55. There is not 
a clear abundance pattern when 2 dimensions are considered ei- 
ther We cannot rule out any projection effects that could blur a 
chemical gradient due to the high inclination of NGC 55. 



5. Summary and Conclusions 

Motivated by the necessity of analyzing large samples of opti- 
cal spectra of massive blue stars in an objective way, even for 
the case of low spectral resolution data, we have undertaken the 
steps to implement a grid-based methodology. We first computed 
an extensive model grid with the model atmosphere/line forma- 
tion code FASTWiND. This new grid was specifically designed for 
the spectral analysis of blue supergiant stars of spectral types 09 
to AO. Secondly, we implemented an algorithm that determines 
the stellar parameters by finding the subset of models in the grid 
fulfilling the criteria of minimizing the differences with respect 
to the observed spectrum. 

We have shown, through a number of control tests, that our 
methodology is well suited for the analysis of optical spectra of 
B-type supergiants, even in the low spectral resolution case. The 
analysis of synthetic models (degraded to the expected obser- 
vational conditions), showed that the stellar parameters are re- 
covered with a high degree of fidelity. Furthermore, the study 



of three Galactic stars, degraded to low resolution and SNR, 
provided answers that are consistent with results based on high 
spectral resolution analysis present in the literature. 

As a first application of our model grid and analysis algo- 
rithm we have analyzed a sample of 12 early B-type supergiants 
located in the Sculptor filament galaxy NGC 55. Our method- 
ology allowed us to obtain a complete characterization of these 
stars, in terms of their stellar parameters and surface chemical 
composition, in spite of the low spectral resolution and SNR. 
The tailored final models provided an accurate match to the ob- 
servations. 

Half of the objects in our sample presented a good agree- 
ment between the evolutionary and spectroscopic masses. For 
the rest, the agreement is not so good, but the results could be 
considered in agreement within the uncertainties of the analysis. 
The location of the stars in the FGLR, for the adopted distance 
to NGC 55, showed a good correspondence with results obtained 
in previous studies. 

The average metallicity of the sample is log (Z/Zq) ~ -0.40 
dex. Our results indicate that NGC 55 does not sustain radial 
abundance gradients, thus confirming previous works based on 
HII regions. Nonetheless, the inclination and morphological 
structure of NGC 55 make this galaxy an interesting target for 
studies of the 2D metallicity distribution, key for understand- 
ing the chemical evolution of NGC 55. The derived CNO com- 
positions show that our stars are evolving away from the Main 
Sequence, and that none of these objects is returning from an 
excursion to the red side of the HR diagram. We have found an 
apparent separation in the nitrogen abundances in two groups, 
both being strongly enriched in comparison with the N baseline 
abundance defined by the HII regions. The derived values are 
in good agreement with a rec ent study of LMC O type stars by 
iRivero Gonzalez et al.l (l20I2h . strongly supporting the idea that 
our objects are evolving directly from the Main Sequence. 

We have shown the reliability of our new automatic, objec- 
tive and fast methodology for the analysis of massive blue stars, 
even at low spectral resolution. Its application to large samples 
will enable us to tackle different issues in a statistical and sys- 
tematic way. In future work, we will apply this method to an 
extended sample of massive blue stars in NGC 55. This analysis 
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Fig. 12. Spatial distribution of chemical species in NGC 55. The plots in the left column display the abundance of Si , Mg , O and 
Z with respect to solar values (Asplund et al. 2009) against projected galactocentric distance {pip,,)- The mean value is indicated 
with a dotted black line and the linear regression with a red dashed line. The linear regressions do not show any clear abundance 
pattern with almost null slopes (see Table |9]l. For oxygen, we also show the H II regions oxygen abundances (blue dots) derived by 
IWebster & SmithI (1 19831) based on the application of strong-line methods. New estimation of O abundances, obtained in this work, 
based on the auroral line [O III] 4363 A measured bv lWebster & SmithI (Il983h (green stars, see text for additional details). The plots 
Ji? the right column show the 2D distribution maps, following the color-code shown in the top-right plot. 
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will provide us with additional information on the discrepancy 
of masses for B-type supergiant stars, their evolution, the FGLR 
and a detailed description of the 2D galactic chemical distribu- 
tion. 
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